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Abstract. We report on a detailed spectral analysis of all the available XMM-Newton 
data of RX J1856. 5-3754, the brightest and most extensively observed nearby, ther- 
mally emitting neutron star. Very small variations (~l-2%) in the single-blackbody 
temperature are detected, but are probably due to an instrumental effect, since they 
correlate with the position of the source on the detector. Restricting the analysis to a 
homogeneous subset of observations, with the source at the same detector position, we 
place strong limits on possible spectral or flux variations from March 2005 to present- 
day. A slightly higher temperature (kT~61.5 eV, compared to the average value kT~61 
eV) was instead measured in April 2002. If this difference is not of instrumental ori- 
gin, it implies a rate of variation of about 0.15 eV yr~' between April 2002 and March 
2005. The high-statistics spectrum from the selected observations is well fit by the sum 
of two blackbody models, which extrapolate to an optical flux level in agreement with 
the observed value. 



1. Introduction 



The X-ray Dim Isolated Neutron Stars (XDINSs, see lTurolla (2009) for a review) are a 



small group of nearby, d < 300 pc, thermally-emitting neutron stars, characterized by 
temperatures kT°° ~ 50- 100 eV, luminosities Lx ~ 10 31 - 10 32 erg s -1 , and spin periods 
in the 3-12 s range. They are radio-quiet and have very faint optical/UV counterparts, 
my ~ 26 - 27. Their spin-down rates (~ 10~ 14 - 10~ 13 ss" 1 ) imply magnetic fields 
B ~ 10 13 - 10 14 G, in good agreement with those inferred from the broad spectral 
features observed in some of these sources. 

Only one of th e XDINSs, RX J0720.4- 3125, showed s ignificant long term spectral 



and flux variations (Ide Vries et al.H 2004; Hoh le et al.l l2012). for which several interpre 



tations were proposed, including precession of the neu tron star and changes induced by 
a glitch (Habe rl et al.ll200d : Ivan Kerkwijk et all 12007 ). No evidence of variability has 



been reported in the other six members of this class, but most of them have not been 
observed very frequently and have lower fluxes than RX J0720.4-3125, therefore the 
derived limits are not very constraining. 
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On the other hand, RX J 1856.5 -3754 is the brightest XDINS and it has been 
routinely observed, almost twice per year since 2002 for calibration purposes, by the 
XMM-Newton X-ray satellite. The resulting large amount of data allows us to inves- 
tigate its spectral evolution on time scales from months to ~10 years, but a reliable 
interpretation of these data requires to carefully consider the stability of the detectors 
and any calibration issue that might affect the results. 
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Figure 1 . Long term evolution of the BB temperature from a set of homogeneous 
observations with the source at the same position on the EPIC pn detector. The dotted 
line is a linear fit to all the points. The dashed line represents a linear fit excluding 
the first point (April 2002). 



2. Results 

We studied the spectral evolution of RX J1856.5-3754 using all the available data ob- 
tained with the pn camera of the EPIC instrument on board XMM-Newton. These con- 
sist of 18 observations carried out between April 2002 and April 2012, for a total expo- 
sure time of more than 700 ks (after screening to remove periods of high background). 
Complex physical models h ave been used in the past t o stud y the X-ray emission from 
RX J1856.5-3754 (see, e.g JTurolla eTaII(l2004l) : lHo1(l2007h ). but, considering the lack 
of prominent spectral features^, and the fact that we are mostly interested in the relative 
spectral variations, we restricted our analysis to fits with one or two absorbed blackbody 
(BB) models. All the fits were performed in the 0.15-1.2 keV energy range. 

The fits with a single BB resulted in temperatures within ~2% of the average 
value kT~61.5 eV. Although the very small statistical errors on the derived tempera- 
tures would imply significant variability, a closer investigation indicate d that the ap 



parent variations are most likely the results of instrumental effects (see iSartore et al. 



'in this respect RX J1856.5-3754 is unique among the XDINSs 



Long term variability of RX J 1856. 5 -37 54 
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Figure 2. Spectrum of RX J 1856.5 -3754 extracted from a homogenous set of 
observations yielding a total exposure time of ~ 254 ks. The upper panel shows 
the best fit with two blackbodies. The features in the residuals (lower panel) have a 
width smaller than the pn energy resolution and are most likely due to non perfect 
instrumental calibrations. 



(|2012l) for details). Due to the varying satellite observing constraints, the source was 
observed at slightly different positions of the EPIC pn detector. We found that the de- 
rived temperature correlates with the source position in instrumental coordinates, in 
particular with the CCD column. This indicates that the effect is probably caused by a 
non-corrected spatial dependence of the channel-to-energy conversion. Note that this is 
a very small effect, which is relevant in this particular case of a high-statistics spectrum 
of a very soft source, and it does not affect most EPIC/XMM-Newton observations. 

To minimize this instrumental effect, we restricted all the further analysis to a 
subset of eight observations in which RX J1856.5-3754 was located at approximately 
the same detector coordinates. In this way we could set strong constraints on the rela- 
tive variations in spectral shape and flux, but it must be remembered that a systematic 
uncertainty of ~l-2% on the absolute values cannot be eliminated. 

The single BB model temperature is compatible with a constant during the last ~5 
years. A linear fit gives a rate of temperature variation of ~ 0.023 ± 0.015 eV yr" 1 , con- 
sistent with a constant at the 2cr level (see FigfT]). A higher temperature was measured 
in April 2002. If not caused by subtle alterations of the instrument response, this differ- 
ence would imply that also RX J1856. 5-3754 undergoes spectral changes, albeit much 
smaller than that observed in RX J0720.4-3125, which underwent substantial and con- 
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tinuous changes in its spectral properties during the years. The temperature variation of 
RX J1856. 5-3754, ~0.5 eV in three years Q corresponds to a rate of - -0.15 eVyr -1 . 

The high-statistics spectrum (~ 2 x 10 6 counts) obtained by summing the data of 
the homogenous observations is shown in Fig. |2j with the two BB best fit model (a 
single BB gives a significantly worse fit). The hard BB has a temperature of kT™ = 
62.4+°-^ eV with emission radius of R™ = 4.7+°j(d/120pc)km, while the soft BB 
has KTf = 38.9^| eV and R™ = 11.8+^ (d/ 120 pc) km. The column density is 
Nu = (12.9 ± 2.2) x 10 1 9 cm" 2 . We found n o convincing evidence for broad or narrow 
absorption features (see lSartore et al.l(l2012h for details). 

The contribution of the soft BB at optical wavelengths is four times larger than that 
of the hard BB. This implies an overall increase of a factor ~5 with respect to the optical 
flux expected from the Rayleigh- Jeans tail of the hard BB alone. The re sulting value 
is con sistent, within the uncertainties, with the observed optical/UV flux (Kapla n et al.1 
20lTh . 



3. Conclusions 

We have carried out a detailed analysis to investigate the long term spectral variability 
of the brightest XDINS, reaching the limit set by the current uncertainties in the in- 
strumental calibration, and without finding evidence for relative variations larger than 
a few percent during the last decade. We stress that the absolute values of the spectral 
parameters derived with XMM-Newton for RX J 1856.5 -3754 (and also for other very 
soft X-ray sources) must be taken with caution since they can be affected by systematics 
errors at the level of few percent. 
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2 or less, considering the lack of coverage between April 2002 and March 2005 



